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This article discusses open-loop frequency-acquisition performance for sup- 
pressed-carrier binary phase shift keyed signals in terms of the probability of de- 
tecting the carrier frequency offset when the arms of the Costas loop detector have 
one-pole filters. The approach, which does not require symbol timing , uses fast 
Fourier transforms (FFTs) to detect the carrier frequency offset. The detection 
probability , which depends on both the 3-d B arm filter bandwidth and the received 
symbol sign ah to- noise ratio, is derived and is shown to be independent of symbol 
timing. It is shown that the performance of this technique is slightly better than 
other open-loop acquisition techniques which use integrators in the arms and whose 
detection performance varies with symbol timing . 


I. Introduction 

Acquiring and tracking binary phase shift keyed 
(BPSK) signals in the absence of a residual carrier is one 
of the many functional requirements that must be met 
by the Block V receiver [1], the Deep Space Network’s 
(DSN’s) next-generation receiver presently under develop- 
ment. The question of how to reduce the frequency error 
A / between the received and predicted carrier frequency 
to within the pull-in range of the Block V receiver’s digi- 
tal Costas loop (the loop used to demodulate BPSK sig- 
nals) was initially addressed in [2]. The open-loop acqui- 
sition techniques discussed in [2] estimate A / by perform- 
ing fast Fourier transforms (FFTs) on the phase detector 
output of Costas-type loops. Because the techniques in 
[2] use integrate-and-dump arm filters, their performance 
in terms of probability of detecting A / depends strongly 
on symbol timing errors, that is, the receiver’s estimate 
of where a symbol epoch starts and ends. In particular, 
the performance of these techniques, which worsens with 


increasing symbol timing errors, motivates evaluating sys- 
tem performance when the integrate-and-dump arm filters 
are replaced by one-pole arm filters that are independent 
of symbol timing. This article considers only the estima- 
tion of the initial frequency error and does not treat the 
acquisition of the phase and frequency after the frequency 
error has been removed and the loop is closed (see Fig. 1). 
Frequency and phase acquisition when the frequency er- 
ror is less than one-half the closed-loop noise bandwidth 
should occur within a few inverse closed-loop noise band- 
widths [3,4]. 

Figure 1 is a functional block diagram of the open-loop 
frequency-acquisition technique under consideration. The 
arm filters are assumed to be one-pole filters. The error 
signal z'(t) is sampled every T sec (T sec is the duration 
of a symbol), accumulated, and fast-Fourier-transformed 
to obtain the error signal spectrum. The probability of 
detecting a tone in white Gaussian noise is well known [5] 
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and can easily be computed as a function of the error- 
signal SNR. The results of [5] can be applied in a straight- 
forward manner because, as shown in Appendix A, the 
error sequence z(n) out of the integrate-and-dump filter 
is composed of a sinusoid at 2A / Hz plus noise which is 
approximately white and approximately Gaussian. The 
absolute frequency difference |2A/| can be estimated (it 
is assumed in the analysis to follow that A fT <C 1) after 
detecting the error signal in the spectrum at the output of 
the Costas loop phase detector. Only the absolute value of 
2A/ can be estimated because taking the product of the 
I and Q arms doubles the frequency error and removes 
its sign. The sign ambiguity problem can be resolved by 
offsetting the frequency so that the error A / always has a 
known sign. 

This article only considers the case of detecting and 
subsequently estimating the frequency error when the er- 
ror signal is restricted to a single FFT bin. It also assumes 
that the error signal is always present. Consequently, the 
article does not address the question of how to choose a 
threshold, which is the power level that the FFT bin with 
the maximum power must exceed before a signal can be 
declared to be present. 

When a tone is known to be present and restricted to 
a single bin, the maximum likelihood (ML) estimate of 
the tone is to observe the FFT magnitude spectrum and 
to select the frequency bin with the maximum power [5]. 
Note that, in practice, it is unlikely that the error signal 
will be restricted to a single FFT bin. 

In a typical deep-space operation, the Doppler shift is 
accompanied by a Doppler rate which is removed by ramp- 
ing the local oscillator phase based on Doppler rate pre- 
dicts. If the Doppler rate cannot be predicted precisely, 
the error signal will drift in frequency during the FFT 
observation time, causing it to smear over multiple bins. 
Clearly, the extent of the smearing depends on the size of 
the Doppler rate error. On the other hand, if the Doppler 
rate is known very accurately so that the signal does not 
drift appreciably during the observation interval, the signal 
may be present in multiple bins due to not sampling the 
discrete-time Fourier transform at its peak (see Chapter 8 
of [6]). Restricting the error signal to a single bin greatly 
simplifies the analysis. Furthermore, it is appropriate be- 
cause the intent is not to quantify absolute performance 
but rather to compare performance between the lowpass 
technique and the integration techniques of [2]. 

Although other models are not precluded, sampling of 
z f (t) is modeled as a T - sec integrate-and-dump filter, be- 
cause doing so defines the signal-tonoise ratio (SNR) of 


the error sequence z(n) in the same bandwidth as the al- 
ternative acquisition techniques studied in [2]. Defining 
SNRs in the same bandwidth is important when compar- 
ing performance between different schemes. Since the arm 
filters do not require symbol timing^ and taking the prod- 
uct of the I and Q arms prior to the integrate-and-dump 
filter gives rise to a signal component (a tone at 2A/) 
whose amplitude does not depend on symbol timing, the 
performance of the one-pole technique is independent of 
symbol timing errors. For convenience, the integrator in 
Fig. 1 is shown as synchronous with symbol transitions al- 
though the analysis does not depend on this fact. In fact, 
the value of integration time T' does not have to be equal 
to T as long as A fV 1. However, as shown in Sec- 

tion II, performance depends directly on the normalized 
bandwidth R , defined as the ratio of the 3-dB arm filter 
bandwidth to the data rate and the symbol SNR. 

A mathematical model of the one-pole technique is de- 
veloped in Section II. Its performance in terms of prob- 
ability of detecting the frequency difference between the 
received and predicted carrier frequency is derived in Sec- 
tion III and discussed in Section IV. Conclusions are stated 
in Section V. 

II. Mathematical Model 

The received sYippressed-carrier BPSK signal, down con- 
verted to an appropriate intermediate frequency (IF), can 
be modeled as 

r(t) = y/2Pd{t)sm(wit + 9i) + n(t) (1) 

where P is the received power in watts, u>i is the IF radian 
frequency expressed in rad /sec, Oi is the signal phase in 
rad, and d(t) is the transmitted data stream given by 

+ oo 

d(t)= d k p(t-kT) (2) 

k — oo 

where p(t) is the baseband non-return-to-zero (NRZ) pulse 
limited to T sec and d k represents the equally likely 
±1 binary symbols. The narrow-band noise process n(t) 
can be expressed as 

n(t) = V2n c {t) cos (ujit + Oi) — \/2n 3 (<) sin(u>,2 -j- #, ) 

( 3 ) 

where n c (t) and n s (t) are statistically independent station- 
ary band-limited white Gaussian noise processes with one- 
sided spectral density No watts/Hz and one-sided band- 
width W Hz. The signal r(t) is demodulated by in- 
phase and quadrature references, a/ 2 sin (uj 0 t + 0 o ) and 
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\/2cos(w 0 i 4- #o), tuned to the predicted carrier frequency 
and then filtered by a lowpass filter. Neglecting higher 
frequency terms, the demodulated signals i(^) and q(t) in 
Fig. 1 can be represented as (see Chapter 5 of [7] for an 
example): 

i(t) = Re ^/Pd(t)e j,p e jAwt ^ - Re n e (t) 

- Re [n J (<)e J ^e- ,Aw( ] (4) 

and 

q(t) = Re ]^/Pd{t) je J ' V Aw ‘ + Re [n c (t)e J ' V Au,t ] 

-Re ^n 3 (t)je>+e> Aut (5) 

where Re is the real part of a complex number, Aw = 
2ttA/ (A/ = fi — f Q ) is the radian frequency error to 
be estimated, and <j> = 0^ — 0 o is the phase error. Let 
£>(w), N 3 (lj) } and N c ( w) be the Fourier transforms of the 
baseband signals d(^), n s (t ), and n c (t). Then the Fourier 
representations of i(t) and q(t) are given by 

I(w) = ^ [o(w - Aw)e j * + ZT (-w - Au)e -J *] 

- ^ [w c (w - A tuy* - N*(-uj - Aw)e - ^] 

- t [Ar 5 (w - Aw)e^ + jV;(-w - Auje'^j 

( 6 ) 

and 

<2(w) = ^ [z)(u> - Aw)e j,A - D*(-w - Auje'^] 

+ ^ [^c(w - Aw)e^ + Af c *(-w - Aw)e‘^l 

- L [at s ( w - A«)e^ - AT;(-w - Aw)e-^| 

(T) 

where A* is the conjugate of the complex number X . Let 
//(w) be the transfer function of the arm filters in Fig. 1. 


Then, the filter outputs in the inphase and quadrature 
arms are given as 

Z,(w) — [-^(w — Ac 4 L>*(— w — Aw)e“ J '*j 

~ — |./V c (w — Aw ) e J * — N c (— w — Aw)e“ J ^J 

- 1 - A«)e^ + 7V/(-w - Aw)c-j* 

( 8 ) 

and 

7,(w) = ^ [l»(w - Au)e^ - D*(-u> - Aw)e‘ J '*] 

+ t j4(w — Aw)e^ + 7V C (— u) — Aw)e“^j 

- Yj [^»( w - - N,’{-u - Aw)e’^] 

(9) 

where D(w), jV c (w), and 7V a (w) are defined, as in [7], as 
D( w) — D(u))H(u) 4 Aw) (10) 

N c {u) = N e (w)H(u>) (11) 

W,(o;) = N a {u>)H{w) (12) 

Equations (8) and (9) represent the output of the arm 

filters when the signal portion of the input has band- 
pass spectra centered at Aw, and the input noise is band- 
limited white noise with bandwidth W B 4 Af where 
B Hz is the 3-dB bandwidth of H(w). Equations (8) and 
(9) are applicable since the demodulated signals ?(<) and 
q(t) have signal spectra centered at Aw during acquisition. 
The time-domain representations of Z,(w) and Z ? (w) are 
given by 

*i(0 = Re [VPd(t)e> V*"*] + ri,(t) (13) 

and 

z q (t) = Re Wpd(t)-. e i*ei^' + n q (t) (14) 
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where 


Ui(t) — —Re 


n c {t)je j *e>* wt 


- Re 


(15) 


n q {t) = Re [n c (ty*e jAwt '. 


— Re 


nAt)\e j *ei Ciwt 

J 


(16) 


and where d(t), n s (t ), and n c (t) are the inverse Fourier 
transforms of D(w), N s {t), and N c {t)- If can be shown 
from Eq. (10) that when Aw ^ 0, d(t) is in general com- 
plex; when Aw = 0, it is real. 


III. One-Pole Arm Filters 

Assuming that the arm filters in Fig. 1 are one-pole 
filters with 3-dB bandwidth B Hz, the transfer function 
//(w) is given by 


H(v) 


U)B 

ju> 4- Wfi 


(17) 


where wb = 2? tB is the 3-dB filter bandwidth in radians. 
For such a filter, the error signal z'(t) — Zi(t)z q (t) can be 
computed as (see Appendix A): 


z'(t) = A(P , AwT, R) sin(2A w* + 2^ + {) + <// (0 

(18) 


and the noise due to the noise times the noise product, 
n f nn (t) } is given by Eq. (A-9). The signal-signal noise, 
n' 4 (<), which is a consequence of intersymbol interference 
(ISI), has two terms. The first term has a continuous spec- 
trum and zero mean, whereas the second term gives rise to 
line spectra at harmonics (not including the fundamental 
harmonic) of the symbol rate. The noise n* SYl {t) is the low- 
pass filtered signal response at time t (due to the present 
as well as previously transmitted pulses) times the thermal 
noise filtered by a lowpass filter. Lastly, the noise-noise 
process n^ n (£) is the product of the filtered thermal noise 
in the inphase arm and the quadrature-phase arm. 

These noises are independent of each other since data 
and noise are assumed to be independent. Consequently, 
the average power of the effective noise is the sum of the 
average power of each of the noise processes above. Ex- 
pressions for average noise power are given by evaluat- 
ing the autocorrelation functions of Eqs. (A-ll) through 
(A-13) at r = 0. 

Referring to Fig. 1, the process z f (t) is integrated and 
dumped to obtain the sequence z(n). Assuming that the 
frequency error is much smaller than the data rate, that 
is, A fT < 1, it is shown in Eq. (A- 15) that the error 
sequence z{n) can be represented as 

z(n) = A(P, R) sin [2Aw(nT + Tf 2) + 2«H- £] + n ef j (n) 

( 20 ) 

where A(P,R) is as defined in Eq. (A-5) and the effective 
noise sequence n e //(n) is defined to be 

««//(") — Uss (n) -j- n an (n) 4- Unn(^) (21) 


The amplitude A given by Eq. (A-5) is a function of many 
variables, including the normalized bandwidth R ~ BT , 
which is the ratio of the 3-dB bandwidth of the one-pole 
filter to the data rate R s = 1/T. It can be shown that A 
does not change significantly with A fT when A fT <C 1. 
Since this article only considers the case A fT <C 1, the 
dependence of A on AwT is omitted in subsequent expres- 
sions to allow a simpler notation. Since the interest here is 
the magnitude of the error signal, the phase £ in Eq. (18) 
is not relevant and consequently is not included here. The 
effective noise into the integrator is defined as 

<//(<) = + + <„(*) (19) 

where the self-noise due to the signal times the signal prod- 
uct, n' 5 (J), ls gi ven by Eq. ( A-7) ; the noise due to the sig- 
nal times the noise product, n' n (i), is given by Eq. (A-8); 


where the noises n 5J (n), n,„(n), and n nn (n) are given by 
Eqs. (A-16) through (A-18). The discrete autocorrelation 
functions of these noises, assuming Aw = 0, are derived 
in Appendix C and listed in Eqs. (A-19) through (A-21). 
These functions are exact except for the autocorrelation 
of n 5a (n), which neglects the negligible second term in 
Eq. (A-7). Since data and noise are independent, the ef- 
fective noise power is given as <r^, where 

4// = k + + *L (22) 

and where Eqs. (A-19) through (A-21) for R.,,(m), 
R sn (m), and R nn (m) yield 

2 P 2 / e - 2 * R + e 2wR \ 2 ( e~^ R \ 

~ 8 ?r 2 P 2 V + 2 ) \ 1 — e -4 *^ ) 

(23) 
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(24) 


<7 


2 

nn 





(25) 


where < 7 2 = N 0 /2T. The signal-to-noise ratio of the error 
sequence z(n), defined as signal power divided by the noise 
power, is given by 


SNR Z (R) = A l( P 2 ’ (26) 

Z<T efJ 

where, after one uses Eq. (A-5) for A(P , /?) and 
Eq. (22) for <r e ^, the term SNR Z can be written in terms 
of the received symbol SNR as 


SNR Z (R) = 


\Co(R)l 

SC ss (R) d- J^TRCsn(R) + JjfftjCnniR) 


(27) 


where SNR = PT/N 0 , |C 0 (i?)| 2 is given by Eq. (A-6), 
and where (7,5(7?) = cr 2 ,/P 2 , C sn (R) = a 2 sn /Pa 2 , and 
C'nn(R) = (Tn n /<7 4 . The quantities C , ,„(/?), and 

C'nr, ( 72) depend on R because cr 2 s , cr 2 n , and a 2 n are func- 
tions of R. 


The objective of the analysis presented here is to enable 
the designer to choose an arm filter bandwidth R for this 
scheme that optimizes the probability of detecting a tone 
in white Gaussian noise, that is, detecting the frequency 
error between the incoming signal and the local oscillator 
frequency so that the reference frequency can be moved 
to the input frequency. Then the loop can be closed to 
speed acquisition. This is accomplished by choosing an R 
that optimizes the SNR of a tone in white Gaussian noise. 
The error sequence SNR given by Eq. (27) represents the 
SNR of a tone imbedded in noise n e //(n), which is neither 
exactly white nor exactly Gaussian but can be assumed to 
be both in practice. This is because the random variable 
in each FFT bin, which results from summing N appropri- 
ately weighted random variables at the FFT input, tends 
toward a Gaussian random variable for large values of N. 
The correlation coefficient p for a one-symbol separation 
is less than 0.2 for symbol SNRs below 0 dB and R > 0.3. 
When p < 0.2, the assumption of independent samples 
out of the integrate-and-dump is valid. Also, it is true 
that n e jj(n) is essentially white when R > 0.5 for SNRs 
above 0 dB. Thus, for any optimum R it is essentially true 
that the integrate-and-dump output sequence is white and 
so the results of [5] apply for the probability of correct de- 
tection of the frequency error. 


IV. Numerical Results and Discussion 

A. SNR Degradation 

Figure 2 depicts SNR degradation D versus normalized 
bandwidth R . Degradation is defined as the reduced error 
signal SNR given by Eq. (26) relative to the SNR of the 
error signal of an “ideal” Costas loop. An ideal Costas 
loop has integrate-and-dump arm filters with (magically) 
perfect symbol timing instead of lowpass filters. In [2] it is 
shown that the error signal SNR for an ideal Costas loop 
detector SNRi is given by 


As expected, for a given SNR, SNR Z is primarily a 
function of R. Equation (27) is exact except for C S$ {R), 
which is approximate because cr 2 s is an approximation: 
see derivation of Eq. (A- 19). Although it is not proven, 
the a 2 s given by Eq. (23) is believed to be a slight upper 
bound to the self-noise power for NRZ pulses. At symbol 
SNRs below 0 dB, SNR Z is a very accurate expression 
for the true error sequence SNR. This is because in this 
region the sum of C sn (R) and C nn (R ), which are scaled by 
the inverse of SNR and are exact, dominates C SS (R). For 
SNRs above 0 dB, SNR 2 is believed to be a slight lower 
bound to the true error sequence SNR, because in this case 
C SS {R) dominates the other two terms in the denominator 
of Eq. (27). 


SNRi = 


SNR 

4 + 2/ SNR 


(28) 


where SNR = PT/No is the received symbol signal-to- 
noise ratio. In mathematical terms, the degradation is 
given as 


D(R) = 


SNRi 

SN R Z (R) 


(29) 


Note that degradation defined in this way is a number 
greater than one, which indicates an actual loss. Fig- 
ures 2(a) through 2(e) show degradation D(R) versus a 
relative 3-dB bandwidth R for various values of SNR. As 
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expected, there is an optimum value for R (i.e., arm filter 
bandwidth) that minimizes degradation and thus maxi- 
mizes error signal SNR. Values of R greater or less than 
this optimum value decrease error signal SNR because they 
increase the noise power faster than the error signal power. 
From Fig. 2 it is clear that lower symbol SNRs result in 
smaller optimum R values. 

The reference SNR, SN Ri given by Eq. (28), is plot- 
ted versus symbol SNR in Fig. 3. Figures 3 and 2 relate 
received symbol SNR to the error signal SNR, SNR Z . 


to a maximum of A/b/at/2. As long as the maximum fre- 
quency error is less than one-half of the loop bandwidth 
(that is, N is large enough), and 

AW < (32) 

2“2 v ' 

or equivalently 

A/b/jv < Bl (33) 


B. Probability of Detection 

The outlier probability (missed detection probability) q , 
defined as the probability that the magnitude of any FFT 
noise-only bin exceeds the magnitude of the signal-plus- 
noise bin, is given by q = 1 — p, where p , the probability 
of detecting the correct frequency offset, is given by (M = 

N/ 2): 


rOQ 

p= 2 M(SNR z )yt 

Jo 


- M(SNR,)(y*+l ) 


x I 0 [2My(SNR z )} [l - y 2 e - M ( SNR '^ M dy 

(30) 


where N is the FFT size (and M is one-half the FFT 
size) and Iq{) is the modified Bessel function of the first 
kind. Thus, this probability for p is given by [5] except 
that the factor of M — 1 is replaced by M in the last 
bracketed term. If M were M — 1, the expressions in [5] 
and Eq. (30) would be identical. The results in [5] are 
slightly different from Eq. (30) because the results in [5] are 
derived for the more general complex FFT case. Figure 4 
depicts q versus SNR where SNR corresponds to the signal- 
to-noise ratio at the FFT input (i.e., SNR Z ). Figure 4 
applies only when the noise prior to the FFT operation is 
white and Gaussian. In the case described here, the error 
signal noise component n e jf(n) is essentially white and 
Gaussian. Since the JV-point FFT has input samples at 
the symbol rate, the frequency bin size is 

A /bin = jjj; (31) 

with (1 jT) as the integrate-and-dump rate which is as- 
sumed to be the symbol rate in this analysis in order to 
compare the results with [2]. Thus, when the correct fre- 
quency bin is detected, the actual frequency error, assum- 
ing that the error signal is at the center of a bin, is reduced 


where B i is the one-sided (closed-loop) bandwidth of the 
Costas loop, the frequency and phase pull-in should be on 
the order of a few inverse Bl s. 

The following example illustrates how to use the curves 
presented in this section to compute the probability of cor- 
rectly detecting the frequency error “seen” by the loop. 
Suppose that the received symbol SNR is 0 dB. Then, 
from Fig. 3, the error signal SNR of an ideal Costas 
loop is —II dB. Assuming that R = 0.3, from Fig. 2(c), 
D( 0.3) = 0.3 dB and the error signal SNR SNR Z = 
-11.3 dB. Finally, from Fig. 4, the probability of incor- 
rectly detecting the actual frequency error for a 1024-point 
FFT is 1.3 x 1(T 6 . 


V. Conclusion 

This article has described a method of determining the 
probability of correctly identifying the frequency error be- 
tween the incoming suppressed-carrier signal and the fre- 
quency of the Costas loop oscillator in order to aid in the 
frequency acquisition of the suppressed-carrier signal. The 
detector chosen for estimating the frequency error is the 
error detector of the Costas loop, which is used for track- 
ing the suppressed-carrier signal. The error signal is not 
fed back to the loop filter and numerically controlled os- 
cillator (NCO), but rather sent to an V-point FFT. The 
FFT then estimates the frequency error, the Costas loop 
is closed, and the loop NCO is adjusted in frequency to 
reduce the initial frequency error to a very small value. 

Knowing the symbol SNR allows the determination of 
the SN R{ and, from Figs. 2(a) through 2(e), the additional 
degradation D. The addition of these two terms yields 
SN R z , which is the abscissa entry on the plot of the prob- 
ability of incorrectly detecting the actual frequency error 
(Fig. 4). One minus this probability yields the probability 
of correctly detecting the initial frequency error between 
the received signal and the rest frequency of the Costas 
loop NCO. 
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This scheme compares favorably to the three methods 
suggested in [2]. In particular, at SNR s that are < 0 dB, 
the staggered approach of [2] is best (minimum degrada- 
tion). Comparing the average loss (averaged over timing 
error) at SNRs < 0 dB, the one-pole arm filter approach 
is about 0.5 dB better. At SNRs > 0 dB it is better by 
more than 0.5 dB. 


Whether the one-pole filter approach or one of the 
methods suggested in [2] will be used for the Block V re- 
ceiver depends on hardware considerations, since the per- 
formance of the best integrate-and-dump arm filter tech- 
nique is nearly comparable to the one-pole arm filter ap- 
proach and is easily switched to the optimum tracking con- 
figuration with integrate-and-dump arm filters. 
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Fig. 1. Costas loop error detector to determine and correct the frequency error between the input and the 

reference oscillators. 
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Appendix A 

Derivation of the Error Sequence z(n) 


Figure 1 shows how lowpass filtered signals are collapsed into an error signal that can be Fourier transformed to 
estimate the frequency error A/. The received signal r(t) downconverted to an appropriate IF is given by Eqs. (l)-(3). 
The signal r(t) is first demodulated and then lowpass filtered to produce the functions Z{(t) and z q (t), which are given 
by Eqs. (13) and (14). Referring to Fig. 1, the error signal z'(t) is given by 

z\t) = Zi(t)z q {t) 


= 1 771 


£<P(t)e , ' a(A " ,+ * ) 


+ Re 


\ypd{ty 




Ug(t) + Re 


\/pJ(O-e ; ' (A “ 1+0) 

3 


n,(t) + hi(t)h q (t) 


(A-l) 


where Re and 7m are, respectively, the real and imaginary parts of a complex number, and where n,{t) and h q (t) are 
given by Eqs. (15) and (16). The quantity d~(t) is given as 

J 2 (0= |d 2 (0|e^ (<) (A-2) 


It can be shown from Eq. (10) that d 2 (i), which is generally a complex quantity, reduces to a real quantity when 
Au = 0. In this case, Eq. (A-l) reduces to (Pd 2 (t)/2) sin 2<f> plus noise. This is the well-known form for the error signal 
of a Costas loop [8], where the quantity d(t) is the filtered version of the baseband data d(t). When Au; ^ 0, the function 
d~(t) for a one-pole filter is derived in Appendix B to be 

OO CO CO 

d 2 (t) = y c n (R)e* 2n * , ' T ) + Y H dkd - lT ^ 1 - 4T ) ( a - 3 ) 

n~ — oo Jb = — co i = -oo 

I** 

where p(t) and C n (R) are respectively given by Eqs. (B-8) and (B-9). Using Eq. (A-3) for d 2 (t) and Eq. (B-4) for d(t) 
in Eq. (A-l), the error signal can be written as 

z f (t) = A(P , II) sin(2A Ljt + 20 + 0 + n' efJ {t) (A-4) 


where 


A(P,R)= -|C 0 («)| 


(A-5) 


id where, from Eq. (B-10), the squared magnitude of Co(R) can be computed to be 


\C 0 (R )\ 2 = 


i + mr 


i+ 


i +(%%y 


(2 TTft) 


7 (1 + e~ 4wR - 2e- 2Tfi cos(A uT)) 


^(l- e -=">(c«s(A„r)-^|sin,i„r) 


(A-6) 
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Since the interest here is in the spectrum of the error signal, the phase £ in Eq. (A-4) is not relevant. Consequently, it 
is not included here. 

The effective noise n' e jj(t) in Eq. (A-4) is defined in Eq. (19), where 


n 'sA l ) = Im 

oo oo 

f E E - kT )p(* - lT)e^ t+ ^ 

+ Im 

p 00 

- C n {R)e R2nTi/T) e j2(Awt+lt,) 

^ n = -oo 


k = — oo i = -°© 
l * k 


n 0 


(A-7) 


”'»«(*) = Re 


VP d k p(t-kT)e j(A “ i+,l>) 


n q (t) + Re 


\/~P dkp(t — kT)—e 
i 3 


j(Au, , f + 0) 


n,(<) 


(A-8) 
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Kni 1 ) = dr{t)n q (t) (A-9) 

where the binary symbols d k are ±1 at random with probability one-half, p(t) and C n (R) are given by Eqs. (B-8) and 
(B-9), and ni(t) and n q (t) are given by Eqs. (15) and (16). These noises are independent of one another since data 
and noise are assumed to be independent. Consequently, the autocorrelation of the effective noise is the sum of the 
autocorrelation of each of the noises. Mathematically, 

K{ f (r) = K a (r) + Kn(r) + R' nn {r) (A-10) 


where, from Appendix C, 
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(All) 


(A- 12) 


(A- 13) 


The preceding equations describe the second-order statistics of the noise, which is usually fed back to close the loop 
along with the signal component . Consequently, the equations can be used to compute the tracking variance of a Costas 
loop with one-pole arm filters. The tracking variance was computed (in the absence of self-noise) by using the noise 
statistics above. It was found to be the same as that computed in Chapter 5 of [7]. Although this is not a definitive test 
of the accuracy of the equations above, it certainly enhances confidence in them. 

As shown in Fig. 1, the error signal z'(£) is integrated and dumped every T sec to obtain the error sequence ~(n). 
That is. 


i An + \ )T 

*(") = / z'm 

J J nT 


(A- 14) 
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Assuming that the frequency error to be estimated, A/, is small relative to the data rate, R s = 1/T, allows one to write 
z(n) as 


z(n) = A(P } R) sin[2A w(nT + Tf 2) + 2<f> + £] + n e// (n) 


(A-15) 


where the sinusoids in Eq. (A-4), which are approximately constant over T sec by assumption, can be removed outside the 
integral after evaluating them at the midpoint of the integration window. The amplitude A(P,R) is given by Eq. (A-5), 
the phase £ is not computed because it is not relevant, and the effective noise n e jj(n) is given by Eq. (21), where 


"««(») = f 

r(n + l)T 
JnT 

(A-16) 

n*n(n) = j; 

,(n + l)T 

/ n ',n(0 rf< 

JnT 

(A- 17) 

n nn {n) = i 

/•{n + l)T 

JnT 

(A-18) 


where and ri f nn (t) are given by Eqs. (A-7) through (A-9). The discrete autocorrelation functions of the 

noises above, which are computed in Appendix C, are 


Rss(m) 


p2 ( , , e~ 2lrR + e 2lrR \ 2 e~^ R 

8^R2V“ 1+ 2 ) I-7--^77. 

0, 


771 — 0 


m/0 


(A- 19) 


Pa 2 


1 L_ ( 1 _ 9 \ , —2ttR ( 1 , 1 _ 17 \ _ e~ 47rJ? , e~ 6?rfi 

1 2vR \2 8irR)J “*■ e y 4 4(2 ttR) S(2tR) 2 J 4(2tt j R) 2 6(2xR) 2 


T 37i (m) = < 


P( 7 2 _L_ fl 2_) + ( 1 + _5_ . . „ 9 \ 1_ (, , _1_\ ( -6* R _ 9,-8 *R , , 

ru 47 tR\ L StvR) ^ 8 ^ 2 ttR^AtFPR 7 ) SttR V 1 ^ SttR M e Ze ^ 


Pcr : 


/ 2 _ e -2^ + 

^ 647r 2 R 2 J \ Z€ 


— 4tt R\m\ _ — 47r/l(|m| + l) _ — 4tt R(\m\— 1 


>). 


m — 0 

10,rfi )] , m = ±l 

m ^ 0,±1 

(A-20) 


^nn(m)= (A-21) 

^ zL ^ e -47r J R(|m|-l) + e -4ir/?(|m| + l) „ 2e“ 4,rH l m l) , TTl ^ 0 

where <r 2 = N Q /2T and i? = BT is the ratio of the arm filter 3-dB bandwidth B to the data rate, 1/T. Note that, 
as R — ► oo, P 33 (rn) — ► 0 as expected. In this case, the data pulses approach the unfiltered case and consequently, the 
self-noise power approaches zero. In the limit R — ► oo, R sn (m) — + Pa 2 , which is the signal-noise power at the output of 
the integrate-and-dump filter when the input is white noise with spectral level PNq/2 watts 2 /Hz. Finally, R n n(m) — ► oo 
as R — ► oo, since in this case the input is white noise with spectral level N^ttB/2, which becomes unbounded as R — ► oo, 
R — BT, for fixed T. 
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Appendix B 

Computing d^(t) in Eq. (A-1) 


The data sequence d(t), which is given mathematically by Eq. (2), has its Fourier transform given by 

oo 

D{u)= d k P(w)e-’ kuiT (B-l) 

k~ — oo 


where the binary symbols dk are ±1 with equal probability and P(u?) is the Fourier transform, denoted ^{}, of the 
baseband NRZ pulse p(t). Applying Eq. (B-l) to Eq. (10) one obtains 

oo 

D(w)= d k P{w)e~> k “ T (B-2) 

Jk = — OO 

where P(w) is defined to be 

P(u>) = P(w)H(u> 4- Aw) (B-3) 

As a result, d(t) = J 7 " 1 {D(u>)} is given as 

oo 

d(t) = d k p(t-kT) (B-4) 

k~ — oo 


where p(t) is the inverse Fourier transform of P(w). From Eq. (B-4), the term d 2 {i) in Eq. (A-1) can be written as 

oo oo oo 

d\t) = J2 p 2 (< - kT) + J2 ^2 d k d,p(t - lT)p(t - kT) (B-5) 


fc— — OO 


k — — oo f = - 00 


The first term in Eq. (B-5) is a periodic complex function with period T. Consequently, it can be represented by a 
Fourier series: 


J2 p 2 {t-kT)= Y, C„e^ 2nrt/T) 

k=-co n = -oo 


where the Fourier coefficients C n are given by 


(B-6) 


Cn = 7pP{p 2 (t)}u,=n(2x/T) 


(B-7) 


For a one-pole filter, p(t) can be computed to be 

p(t) = ^ ( p(t) - e~( WB+j * u)t u(t) + - T)) (B-8) 

u> B + j Aw V ' 

where u> B is the 3-dB filter bandwidth in radians, p(t) — u(t) — u(t — T), and u(<) is the unit step function. Applying 
Eq. (B-8) to Eq. (B-7), the Fourier coefficients C n as a function of A ujT and R — BT are computed to be 
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C n (R) = 




sin nn 
nn 


1 (1 + 0 - e-( 2 '«+^" T )) 

2^(l+j2-^)(i + i -^21) 


(B-9) 


Note that for the special cases of Aw = 0 and A u> = (— n7r)/T, Eq. (B-9) reduces to Table 2 of [9]. Additionally, for 
n = 0, Eq. (B-9) reduces to 


C 0 (R) = 
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Appendix C 

Computing R’ss (t), R’sn (t), R’nn (t) and Rss {m)> Rsn (m), Rnn {m) 

When applicable, the derivations that follow assume that the frequency error A / = 0. This greatly simplifies the 
computation of the various autocorrelation functions and is not expected to change appreciably the final result of the 
main text. The reason for this is that all the significant noise components give rise to continuous spectra whose shapes 
are not expected to change much with A / when A / <C 1/T. Furthermore, when computing fi 33 (r), the second term in 
Eq. (A-7) is not included. This approximation is made because, except for the first harmonic which has a line component 
at 1/T ± 2A /, the remaining harmonics give rise to line components far beyond the band of frequencies within which 
the error signal is expected. Additionally, power in the first harmonic which is at least 8 dB down from the component 
at 2A / when A uiT < 0.1, is suppressed even further by the integrate-and-dump filter response preceding the FFT 
operation. 


I. R'ss (t) and Rss (m) 

The self-noise n 33 (<) is given by Eq. (A-7). When Aw = 0, p(t) is real, and n 33 (t) reduces to 


,«) = f 


E E d k d,p(t - lT)p(t - kT) sin 20 


V 


k OO 00 


(C-l) 


where d n is ±1 at random with probability one-half, and 0, which is independent of d n , is a uniform random variable in 
the interval [0, 2 tt]. The filtered pulse p(t) is given by Eq. (B-8). The autocorrelation function of n' 33 (<) is given by 


K,( r ) = £K SJ (*K*( f + r )] = ~g’ R ( r ) 


(C-2) 


where, after noting that £ [sin 2 <fi] = ^ , R(t) can be defined as follows: 

oo oo oc oo 


R(r)= E E E E d kdld m d n p(t - kT)p(t - lT)p(t - mT)p(t - nT) 


(C-3) 


k = — co l~-°° m = — oo n=-. 

l*k njtr 


The function above, which is in general difficult to compute because the filtered pulse p(t) is not time limited, can be 
approximated as follows. One begins by expressing R(r ) as the inverse Fourier transform of its corresponding spectral 
density. 


R(t) = 



did 

5(w)e-- 


« 5(0) lim 

a;— ►oo 



0 3 ^t 
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From Chapter 5 of [6], the dc component of the spectral density in Eq. (C-4) is given by 


5(0) = 4T R? d (lT) 
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where S d ( u>) is the spectral density of the random NRZ data and |//(w)| 2 is the squared magnitude of the one-pole 
filter with 3-dB bandwidth B IIz. Grouping the exponential in the above integral with \H(u>)\ 2 and applying Parseval’s 
theorem gives the corresponding integral in the transformed domain. Namely, 


Rt[,T) = L (' - r) w (- 


1 + 


e -2*R + e 2*R 


(C-7) 


where the parameter R = BT. Finally, R',,(t ) ' n Eq. (A-ll) follows after back-substituting the last result. 

The discrete autocorrelation function given by Eq. (A-19), which is the correlation function of the noise sequence 
n ss (n) given by Eq. (A-16), follows directly upon evaluating the next equation: 

= —7 / / R'sAlh ~ ti\)dhdti (C-8) 

1 “ J 0 JmT 


where R , ss (r) is given by Eq. (A-ll). 


II. R' sn (t) and R sn ( m ) 

The signal-noise process n' n ($) is given by Eq. (A-8). For A u) = 0, p(t) is real and n' n (<) reduces to 


= f \fP ^ d k p(t - kT) cos <j) h q (t) + >/P ^2 d k p(i — kT) sin </> j n,*(<) 


Jfe=- 


k = - 


(C-9) 


Since the data and noise processes are independent of each other, the autocorrelation of the signal-noise product is the 
product of the autocorrelation of the individual signal and noise processes. Hence, 

KJr) = £[n' an (t)n' in (t + r)] = Pfi J (r)/i n (r) (C-10) 


where 


R,(t) = £ 


OO OO 

E E d k d,p(t - kT)p(t - IT) 
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(C-ll) 


Rn(r) = £[ni{t)hi{i + r)] = £[h t (t)h q (t + r)] (C-12) 

The functions R s {r) and Rn(r) are the autocorrelation functions of signal and noise processes at the output of a one- 

pole fdter when the input signal process is random NRZ data, and the input noise process is white noise. Consequently, 

writing R s (r) as the inverse Fourier transform of its corresponding spectral density in terms of the input spectral density 
times the squared magnitude of the one-pole filter, grouping the exponential with the squared magnitude of the filter 
transfer function, and applying ParsevaFs theorem yields 

R,(r) = j\ (1 - !£!) (C-13) 

where B Hz is the 3-dB bandwidth of the one-pole filter. The noise correlation is well known to be 

R„(r) = (C-14) 
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Finally, the signal-noise correlation given by Eq. (A-12) follows after substituting Eqs. (C-13) and (C-14) into Eq. (C-10). 

The discrete autocorrelation function given by Eq. (A-20), which is the correlation function of the noise sequence 
rc 3 „(n), is obtained upon evaluating the next equation: 

i r T 

Rsn(m) - — / Kn(\t* - h\)dt 2 dti (C-15) 

1 ‘ J0 JmT 


Since the integrand above depends only on the absolute difference 1 1 2 -*i|, the double integral above can be reduced to 
the following simpler single integrals when m = 0, 1, — 1. Namely, 


T f (m+l)T 
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(C- 16) 


When |m| > 1, Eq. (C-15) reduces to the product of two single integrals in a straightforward manner. 


III. R'rw (r) and R nn (m) 

The noise process n' nn (t) given by Eq. (A-9) is the product of the filtered inphase and quadraphase noise processes, 
hi(t) and h q (t ). The autocorrelation of these noises, denoted by # n (r), is given by Eq. (C-14). Furthermore, since 
hi(t) and n q (t) are independent with respect to each other, R nn (r) = Rn( r ) an< ^ Eq. (A- 13) follows. The discrete 
autocorrelation Rnn(m) is given upon evaluating the following equation: 

pT A(m + 1)T 

R U(m)= / / R' nn {\t 2 - U\)dt 2 dt x (C-17) 

Jo J mT 


where, when m ~ 0, the double integral above can be transformed to the following simpler single integral: 



Rnn (1*2 — ti\)dt2(lti 



(M) ( T - M) Jr 


(C- 18) 


When m ^ 0, Eq. (C-17) reduces to the product of two single integrals. 
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